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SUMMARY 


This  report  is  concerned  with  the  mathematical  features  of  the  de- 
sign of  a control  system  that  has  the  property  of  simultaneously  reducing 
the  error  and  the  derivatives  of  the  error  to  zero  in  a minimum  time.  The 
differential  equation  of  the  control  system  is  assumed  to  be  of  the  form 


Ly  = fjy. 


dt  5 ' * ' 


i 


9 


where  L is  a linear  ordinary  differential  operator  of  the  n^  order 
with  constant  coefficients,  y is  the  erro**.,  and  f is  the  control  func- 
tion or  correcting  force.  In  most  of  the  report,  it  is  assumed  that  the 
function  f is  discontinuous  and  takes  on  only  the  values  -1  . Thus, 
the  control  function  acts  as  a simple  switch  or  relay.  The  problem  is  to 
determine  how  the  function  f should  defined,  i.e.,  to  determine  how 
to  switch  the  correcting  force  from  +1  to  -1  and  back,  so  that  the 
error  and  its  derivatives  are  reduced  to  zero  in  a minimum  time. 


The  n -order  equation  is  reduced  to  a system  of  first-order  equa- 
tions and  the  techniques  of  matrices  are  used.  The  problem  is  solved  for 
all  cases  when  she  eigenvalues  of  the  matrix  of  the  system  (or  the  roots 
of  the  characteristic  equation  of  the  operator  L)  are  real.  The  result 
can  be  stated  as  follows:  If  it  is  possible  to  find  any  control  function 

f that  simultaneously  reduces  the  error  and  its  derivatives  to  zero,  then 
there  exists  a unique  function  f which  does  this  in  a minimum  time.  The 
minimizing  function  is  the  one  which  employs  a minimum  number  of  "switches." 
It  should  be  pointed  cut  that  this  result  holds  only  for  real  eigenvalues. 


A slightly  broader  problem  is  considered  in  a few  simple  special 
cases,  namely,  the  problem  where  f is  restricted  only  by  the  inequality 
If  | fC  1.  It  is  shown  in  these  cases  that  no  reduction  in  the  minimum 
time  is  possible,  and  that  the  minimizing  function  is  still  a function  of 
the  discontinuous  type,  where  f = -1  . 


The  study  was  conducted  at  the  Experimental  Towing  Tank,  Stevens 
Institute  of  Technology,  under  Office  of  Naval  Research  Contract  No. 
Nonr-26302. 
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INTRODUCTION 

Limit  control  refers  to  a discontinuous  control  that  is  allowed  to 
take  on  onlv  extreme  values.  Such  controls  are  often  referred  to  as 
bang-bang,  relay  type,  or  on-off  controls.  The  present  study  is  con- 
cerned with  the  design  of  a limit  control  system  that  will  simultaneously  ^ 

reduce  the  error  and  the  derivatives  of  the  error  to  zero  in  a minimum  ! 

time. 

Mathematically,  the  error  is  represented  by  a differential  equation  I 

or  a system  of  differential  equations  where  the  forcing  function  or  con- 
trol function  is  discontinuous  and  is  allowed  to  take  on  only  extreme  | 

values.  The  problem  is  to  find  how  to  switch  the  control  function  from 
one  extreme  value  to  the  other  in  such  a way  that  the  error  and  its  de-  | 

rivatives  are  reduced  to  zero  in  a minimum  time. 

This  report  extends  some  of  the  results  presented  in  Reference  1,  | 

I 

which  were  essentially  the  first  results  for  this  problem.  The  work  is 

part  of  the  limit  control  project  of  the  Experimental  Towing  Tank,  Stevens  | 

Institute  of  Technology,  under  contract  with  the  Office  of  Naval  Research. 

I 

I 

I 
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STATEMENT  OF  THE  FROBLEM 


An  Automatic  Control  Problem 

This  study  is  concerned  with  a minimum  problem  associated  with  a 
linear  differential  equation  containing  a discontinuous  forcing  term.  The 
problem  arises  physically  in  the  design  of  an  automatic  control  system  for 
ontimum  performance. 

Consider  the  control  system  shown  in  Figure  1: 


1 

Controlling  . 

£ 1 

Control  led 

I y 

1 * 

j 

System  j 

1 i 

System 

Power 

Source 


FIGURE  1 

The  purpose  of  this  system  is  to  keep  the  output  y of  the  controlled 
system  constant  at  the  value  y . If,  at  any  instant,  the  output  is 
different  from  y , then  an  error-sensitive  device  feeds  the  error  y - y 
back  to  the  controlling  system.  The  output  of  the  controlling  system,  i.e. 
the  control  function  f , is  fed  into  the  controlled  system  and  should  be 
designed  so  as  to  cause  a reduction  in  the  error.  The  important  part  of 
the  system  is  the  control  function  which  is  supposed  to  be  some  function 
of  the  error  and  its  time  derivatives  up  to  a certain  order.  The  broad 
problem  considered  here  can  now  be  stated: 

To  determine  the  control  function  so  that  the  error  and  the  de- 
rivatives of  the  error  up  to  a certain  order  are  simultaneously  reduced 
to  zero  in  a minimum  time. 

In  what  follows,  a drastic  restriction  is  made  on  the  class  of 
control  functions  that  are  considered.  In  fact,  in  all  but  one  case. 


*-!i59 
- h - 


it  wil]  be  assumed  that  the  control  function  f jL£  a discontinuous  function 
that  is  allowed  to  take  on  only  the  values  ♦b  and  -b  , where  b is  a 
positive  constant/.  The  controlling  system  therefore  acts  as  a switch  or 
a relay  that  applies  the  full  strength  of  the  power  source  to  the  controlled 
system.,  either  directly,  or  with  a reverse  in  polarity. 

There  is  no  conclusive  reason  why  a relay  type  control  should  be  the 

best  one  for  the  present  problem.  However,  since  the  full  strength  of  the 
pow-er  source  is  always  used,  it  does  seem  reasonable  that  the  error  and 
its  derivatives  can  be  reduced  to  zero  in  a shorter  time  than  by  any  other 
means.  In  this  connection,  it  will  later  be  proved,  only  in  the  case  of 
very  simple  systems,  that,  among  all  functions  satisfying  |f  |<  b,  the 
minimizing  function  is  of  the  discontinuous  type  described. 

Mathematical  Formulation  of  the  Problem 

In  the  system  shown  in  Figure  1,  it  is  assumed  for  convenience  that 
the  reference  value  v^  is  zero.  The  error  is  then  simply  y(t)  „ It 
is  also  assumed  that  the  controlled  system  is  governed  by  a linear  dif- 
ferential equation  with  constant  coefficients  so  that  the  error  y(t)  is 
a solution  of 


j ,n-l  \ 

Lv  = r|v  & d jr 

^ dt  n-l  j 3 


where  L represents  the  operator 


LEDn  + a^11'1  + 


. . . a ,D  + a 

n-l  n 


(1) 


(2) 


and  tne  control  function  f is  a discontinuous  function  taking  on  only 
the  values  +1  and  -1  . (if  f = -b  , then  b can  he  made  to  equal  1 
by  a change  in  the  time  scale.)  A preliminary  statement  of  the  problem 
is  to  determine  the  function  f (that  is,  to  determine  for  what  values 
f = +1  and  for  what  values  f ~ -l)  such  that  the  solution  of  equation  (l) 
which  satisfies  the  initial  conditions 


y(0) 


dy(0) 

dt 


^(O)  = 


dt 


n-l 


(c^  - constants). 


will  reach  the  state 


y 


= = 
dt 


= 0 


in  a minimum  time. 

If  the  function  f is  at  all  complicated,  there  is  some  doubt  as  to 
what  is  meant  by  a solution  of  equation  (l)  (suppose,  for  instance,  that 
f is  +1  when  y is  an  irrational  number  and  f is  -1  otherwise). 

To  circumvent  this  difficulty,  the  problem  is  formulated  in  different  terms. 
For  this  purpose,  equation  (l)  is  reduced  to  a system  of  first- order 
equations.  In  the  usual  manner,  let 


t r = 


> (3) 


n-l  n 

= -a  x..  ••  a ,x,.  - 
n n 1 n-l  2 


, . - anx  + f 
1 n 


J 


(x  = dXj/dt  etc.).'  In  matrix  form,  this  can  be  written  as 

x = Ax  + p , (U) 

where 


"xl~ 

‘ 0 

.1 

0 . 

. 0 

0 1 

0 

x2 

0 

0 

1 . 

. 0 

0 

0 

• 

, A = 

• 

* 
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• 

• 

• 

. 

• 

; o 

0 

0 . 

. 0 

1 

0 

X 

n 

i 

-a 

L n 

~an-l 

4 • 

. -a2 

-al_ 

f(x) 

Since  f can  take  on  only  the  values  -1  , the  solutions  of 
equation  (h)  shoula  be  combinations  of  the  solutions  of 
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x = Ax  - e , (5) 

where  e = col(0,0, . . . ,0,l)  . The  solutions  of  x Ax  + e will  be 
called  P-curves  (P  for  positive)  and  those  of  x = Ax  - e will  be  called 
N-curves  (N  for  negative).  Along  each  P-  or  N-curve  there  is  a natural 
orientation,  namely,  the  direction  in  which  t increases. 

/ \ \ o 

A solution  of  (4)  from  the  point  x is  defined  as  a continuous 
curve  that  consists  of  a finite  sequence  of  alternating  P-  and  N-arcs 
(each  of  finite  time  length)  starting  at  the  point  x°  . The  solution 
can  be  represented  by  x(t)  , a continuous  vector  function  of  the  scalar 
parameter  t . A path  from  x°  is  defined  as  a solution  from  x°  which 
ends  at  the  origin.  With  each  path  there  is  associated  a transit  time, 
namely,  the  sum  of  the  time  lengths  of  the  F-  and  N-arcs  making  up  the  path. 

With  this  terminology,  the  problem  can  finally  be  stated  as  follows: 

To  find _a  path  of  minimum  transit  time  (a  minimal  path)  from  each 
point  x°  in  the  phase  space. 

Corresponding  to  each  solution  of  equation  (U),  as  defined  above, 
there  is  a control  function  f , namely,  f = +1  on  the  P-arcs  and  f = -1 
on  the  N-arcs  (here,  a P-  or  an  N-arc  is  regarded  as  being  closed  at  its 
initial  point  and  open  at  its  terminal  point).  The  only  functions  f 
admitted  are  those  such  that  all  the  solutions  from  a point  are  of  the 
type  described. 

Since  the  problem  has  been  formulated  entirely  in  terms  of  the  system 
of  equations,  it  is  not  necessary  to  restrict  the  matrix  A to  be  of  the 
above  form.  The  problem  can-  be  generalized  slightly  and  instead  of  con- 
sidering equation  ($),  the  following  equations  will  be  considered: 

x - Ax  - e , (6) 

where  A is  any  n by  n matrix  and  e is  any  non-null  vector. 


* Superscripts  on  vectors  are  indices  and  not  exponents. 

**  It  shoul^i  be  remarked  that  equation  (6)  is  not  necessarily  equivalent 
to  an  n^- -order  equation.  However,  physical  systems  arise  directly  in 
the  form  (6),  such  as  the  linearized  equations  that  govern  the  motion 
of  a submerged  submarine. 


It- 


[ ' 


i 
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The  final  formulation  of  the  problem  to  be  treated  in  this  report  is: 

Consider  the  system  of  differential  equations 

- — — ^ — — — — ' •' 

x = Ax  + ef(x)  , (7) 

where  A is  any  n by  n matrix,  e is  any  non-null  vector,  and  f(x)  , 
a.  real- valued  function  of  the  vector  x , can  take  on  only  the  values  -1  . 
How  should  the  function  f(x)  be  chosen  so  that  the  solution  x(t)  satisfy- 
ing arbitrary  initial  conditions  x(0)  = x°  reaches  the  origin  x = 0 in 
the  least  possible  time? 


General  Remarks  About  the  P-  and  N-Curves. 

The  existence  or  nonexistence  of  minimal  paths  of  course  depends  on 
the  behavior  of  the  P-  and  N-curves.  Some  of  the  properties  of  these  curves 
are  discussed  below. 

An  important  fact  is  that  the  family  of  P-curves  is  symmetric  to  the 
family  of  H -curves  with  respect  to  the  origin.  This  holds  for  orientation 
as  well  as  for  shape.  The  symmetry  follows  from  the  fact  that  if  x(t) 
satisfies  x = Ax  + e with  x(0)  = x°  , then  z(t)  = -x(t)  satisfies 
z = Az  - e with  z(0)  = -x°  . Using  this  property,  it  is  necessary 
to  solve  the  minimal  problem  for  only  one-half  of  the  phase  space,  the 
solutions  for  the  other  half  being  obtained  by  symmetry. 

It  is  known  from  the  well-known  existence  theorem  for  linear  dif- 
ferential equations  that  there  exists  a unique  P-  and  N-curve  through  each 

point  of  the  phase  space.  The  equations  of  these  curves  can  be  found  in 

A ^ 

terms  of  the  matrix  function  £ . The  solutions  of  the  differential 

equations 

x = Ax  - e (8) 


satisfying  the  initial  conditions  x(0)  = x°  are 

,t 


. At  o + At 


.-As 


x = € x - 6 J € ds  e , 
0 ' 


(9) 


i ' 
! 


which  are  the  equations  of  the  P-  and  N-curves.  (In  the  case  of  an 
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L 


I 
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ambiguous  sign,  the  upper  sign  holds  for  the  P-curves,  the  lower  sign  for 
the  N-curveSo) 

If  the  matrix  A is  ncnsingular,  (9)  can  be  simplified  to 


x 


At  0 + 
= € x - 


I)A“1e 


(10) 


which  can  be  verified  by  direct  substitution.  The  points  +A  e are 
singular  points  (or  equilibrium  points)  of  the  DE  , that  is,  points 
where  x = 0 . The  P-curve  through  the  point  -A  is  just  the  point 
itself  and  similarly  for  the  d -curve  through  +A  e » 

A.'tt 

It  is  well  known  that  the  matrix  £ “ can  be  determined  explicitly 

if  the  eigenvalues  of  the  matrix  A are  kncwn.  In  fact,  in  order  to  solve 

the  problem,  it  will  be  convenient  to  make  a substitution  in  the  DE  that 

will  reduce  the  matrix  A to  diagonal  or  Jordan  canonical  form  where  the 

At 

eigenvalues  are  explicitly  displayed  and  where  the  function  € assumes 
a simple  form.  Obviously,  the  character  of  the  P-  and  N-eurves  is  also 
determined  by  the  eigenvalues  of  A . It  will  be  shown,  in  fact,  that 
the  character  and  existence  of  minimal  paths  are  determined  solely  by  the 
nature  of  the  eigenvalues. 


In  this  report,  a complete  solution  to  the  problan  will  be  given  in 
the  case  where  the  eigenvalues  of  A are  all  real.  The  trivial  one- 
dimensional case  will  be  solved  first;  then  the  two-dimensional  cases 
will  be  discussed,  and  these  will  be  used  as  a basis  for  induction  in 
solving  the  n-dimensional  cases.  In  addition,  the  broader  problem  where 
the  control  function  f is  restricted  only  by  the  inequality  | f | < 1 
will  be  solved  for  the  cases  of  the  one-dimensional  systans  and  a simple 
two-dimensional  system. 


The  only  known  previous  results  for  the  present  problem  are  given 
by  Bushaw  in  Reference  1.  Bus haw  completely  solved  the  problem  for  all 
two-dimensional  systems  (including  the  case  of  complex  eigenvalues)  that 
arise  from  a second  order  differential  equation  with  constant  coefficients. 


Before  proceeding  with  the  analysis,  a few  remarks  will  be  made  about 
the  general  procedure.  The  problem  is  to  find  the  minimal  paths  which,  by 
definition,  must  enter  the  origin.  The  only  way  a minimal  path  can  reach 


the  origin  is  by  way  of  the  P--  or  d -curve  through  the  origin.  Let  the 

curve  T be  defined  as  the  curve  obtained  by  following  the  P-curve  through 
the  origin  backward  in  time  (starting  at  the  origin).  The  minimal  path  for 
at  least  some  of  the  points  of  T close  to  the  origin  must  consist  of  just 
that  portion  of  T between  the  point  and  the  origin  (if  a minimal  path 
exists  at  all).  By  symmetry,  corresponding  points  are  found  on  the  N-curve 
through  the  origin,  Thus,  all  the  minimal  paths  with  no  "corners"  are 
found,  a " corner 11  being  a junction  of  nonzero  P-  and  N-arcs.  Then  the 
minimal  paths  with  exactly  one  corner  are  determined,  and  so  on,  until 
all  che  minimal  paths  have  been  found. 
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ONE-  AND  TWO-DIMENSIONAL  SYSTEMS 


One-Dimensional  System 

The  differential  equations  for  a one-dimen3ional  system  are  simply 
x.^  = ax.^  - « (ll) 


■where  x^  , a , e^  , are  all  real  numbers.  It  is  assumed  that  e^>  0 . 

Three  cases  need  to  be  considered,  namely  a < 0 , a > 0 , and 
a = 0 (a  is  actually  the  eigenvalue).  If  a ^ 0 , then  the  solutions 
of  ( ll)  are 


at/  c + “IN 

e 1*1  ‘T  / 


(12) 


which  are  the  equations  of  the  P-  and  N-curves  whose  initial  point  is  x^0  . 

Consider  first  the  case  a < 0 The  P-curves  start  at  the  point 
x^°  and  continuously  approach  the  point  -e^/a  as  t-»-00  . The  N-curves 
continuously  approach  the  point  +e-^/a  as  t 00  . This  behavior  is 
shown  schematically  in  the  following  Figure: 


FIGURE  2 

How  to  find  the  minimal  paths.  If  x^0  is  in  the  interval 
0 x^  — e.^/ ^ > the  minimax  pain  is  easxxy  seen  uo  be  the  N— arc  from 

x^°  to  the  origin.  This  follows  from  the  fact  that  the  N-arcs  in  this 
interval  are  directed  toward  the  left  while  the  P-arcs  are  directed  toward 
the  right.  If  the  minimal  path  contained  a P-arc,  this  would  mean  that  x^ 
would  be  increased  and,  in  order  to  reach  the  origin,  a longer  N-arc 
would  have  to  be  traversed. 

If  x^°  is  in  the  interval  x^°  > -e,/a  , both  the  N-arcs  and  the 


P-arcs  are  directed  tcward  the  left.  However,  the  minimal  path  is  still 
the  N-arc  connecting  the  point  x, 0 with  the  origin.  Because  of  the  re- 
sult in  the  preceding  paragraph,  it  is  only  neci  sary  to  show  that  the 

minimal  rath  cannot  contain  a P-arc  in  the  interval  x,  > -e.,/a  . If  the 

■l  £ 

minimal  path  did  contain  a P-arc,  say  between  the  points  x^  and  x^  , 

where  xn°  and  x are  greater  than  -e  /a  and  x,  < x 0 , the  time  to 
m X X X X 

traverse  this  arc  could  be  easily  obtained  from  (12)  to  be 


The  time  to  traverse  an  N-are  connecting  the  same  points  is 


(13) 


(lli) 


It  is  easy  to  prove  that  t^  < tp  and  therefore  that  the  minimal  path 
cannot  contain  a P-arc 

By  symmetry,  an  analagous  result  is  obtained  fcr  the  points  x^°  < 0 
where  N-arcs  are  replaced  by  P-arcs.  hence,  a unique  minimal  path  exists 
from  every  initial  point. 

If  the  eigenvalue  a is  positive,  the  equations  for  the  P-  and 
N-curves  are  still  given  by  (li)  but  the  directions  are  reversed,  as  shown 
in  the  following  Figure: 


( 


0 


P 


► 


FIGURE  3 


It  is  seen  immediately  that  no  oaths  exast  if  | x^°  | > e^/a  . If 
| x^-"  j < e^/a  , the  same  procedure  as  above  proves  that  a unique  minimal 


path  exists.  If  0<  x^°<C  e^/s  > the  minimal  path  is  the  N-arc  connecting 
x 0 with  the  origin,  and  if  -e^/a  — xi°  ^ 0 , the  minimal  path  is  the 
F-arc  connecting  x^°  with  the  origin. 

If  a = 0,  the  solution  (1?)  fails.  The  equations  of  the  P-  and 
N-curves  in  this  case  are 


x 


1 


(15) 


Here,  the  P-curve  from  any  point  is  directed  to  the  right  and  the  N-curve 
is  directed  to  the  left.  By  using  the  same  reasoning  as  above,  it  is  seen 
that  a unique  minimal  path  exists  from  every  point  x^°  ; if  >:^°  > 0 , the 
minimal  path  is  the  N-arc  through  x-^0  , while  if  x^°  < 0 , the  minimal 
path  is  the  P-arc  through  x^°  . 

The  results  of  this  section  are  summarized  in  the  following  theorem: 
Theorem  1:  In  the  one-dimensional  system 


x = ax  + o ffr  ) 

i i r ' r 


where  e1  > 0 


and 


f(xx) 


is  restricted  to  the  values 


, a unique 


minimal  path  exists  if  a_<0;if  a > 0 , a unique  minimal  path  exists 
for  points  x^0  such  that  j x^°|  <e^/a  and  no  path  exists  for  | x^°|  > e^/a 
The  minimal  path  (if  it  exists)  for  x^  > 0 is  just  the  N-arc  connecting 
x^°  with  the  origin,  and  for  x-^0  <0  , the  P-arc  connecting  x^°  with 
th3  origin,  (in  terms  of  the  control  function  f , this  means  f(x^)  = 

-s  ignum  x^  . ) 


One-Dimensional  System  with  |f(x^)j  O. 

The  broader  problem  will  be  considered  here  for  the  one-dimensional 

system 

= ax^  + eLf(x  ) , (16) 

where  e^  > 0 and  |f(x^)|  <1.  In  order  to  insure  solutions  of  the 
differential  equation,  it  will  be  assumed  that  f(x^)  is  piecewise 
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continuous  with  only  a finite  number  of  discontinuities  in  any  bounded  in- 
terval. The  solution  of  ( 16)  through  any  initial  point  will  then  be  uniquely 
determined  by  the  sole  requirement  of  continuity.  The  class  of  functions 
considered  here  will  be  broader  than  the  class  of  functions  considered  in 
the  preceding  section.  It  will  be  shown,  however,  that  the  minimum  time 
obtained”  in  the  last  section  cannot  be  improved. 

Theo  rem  2_:  Among  all  functions  f (x^)  satisfying  the  conditions 

in  the  preceding  paragraph,  the  minimizing  function  is  just  the  function 
f(x^)  = -signum  x^  that  was  obtained  in  Theorem  1. 

This  theorem  will  be  proved  only  for  the  case  a < 0 . The  other 
cases  can  be  treated  in  a similar  manner.  If  the  initial  point  is  x^°  >0, 
the  time  T for  fCx.^)  = -signum  x^  to  get  to  the  origin  is  simply 


T = 


/a 


(17) 


It  will  be  shown  that  the  time  for  any  other  f will  exceed  T . Since 
j f | < 1 and  e^>  0 , 

~el  <e1f(x1)  < e±  (18) 


or,  using  the  differential  equation  (16), 

-e1<  xx  - ax1  <ex  • (19) 

-a.tr 

Multiplying  (19)  bv  the  integrating  factor  £ gives 


,-at 

e„c 

x 


Integrating  tnis  inequality  from  0 to  t , 
gives 


(20) 

remembering  that  x^(0)=x^°. 


a 


(21) 


i 
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Setting  x^  = 0 and  solving  the  left-hand  inequality  results  in  the  follow- 
ing inequality  for  the  time  to  reach  the  origin: 


1 X1  " el/a 

-.A 


This  shows  that  t > T . It  is  easy  to  see  that  the  strict  in- 
equality Kill  hold  unless  the  function  f(x^)  is  identical  with  the 
minimizing  function  f(x^)  = -signum  x^  . This  completes  the  proof. 

Two-Dimensional  Systems 

The  minimal  problem  will  be  considered  for  a two-dimensional  system 


x = Ax  - e 


(e  / 0)  , 


where  A is  a 2 by  2 matrix  and  x and  e are  two-dimensional  vectors. 
Onlv  those  cases  where  A has  real  eigenvalues  will  be  considered. 

It  is  convenient  to  make  a substitution  in  the  DE  , namely, 
x = Tz  , (2k) 


where  T is  a real  nonsingular  matrix.  The  DE  becomes 


z = Dz  - w , (25) 

where  D = T "''AT  and  w = T "'"e  . The  matrix  T can  be  chosen  so  that  D 

it 

assumes  one  of  the  three  following  forms: 


Case  (a) : D = 


Case  (b): 


\1  0 
.0  X 


r\ 


0 x 


X 1 


Case  (c):  D = , 

0 X * 


* See,  e.g. , Reference  2,  page  20tf« 
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where  X^  , X are  the  eigenvalues  of  A . Case  (a)  will  be  obtained  if 
the  eigenvalues  are  distinct,  case  (b),  if  the  eigenvalues  are  equal  but 
there  are  two  linearly  independent  eigenvectors,  and  case  (c),  if  the  eigen- 
values are  equal  with  only  one  independent  eigenvector. 


The  solutions  of  (25)  are 


,Dt  o ♦ Dt 
z = € zu  - € 


1 


where,  corresponding  to  each  of  the  forms  of  D in  (26), 


(27) 


Case  (a):  fDt  =,  € 


0 


xt 


.xt 


Case  (b) : 6 


Dt 


0 


„ X t ! 
O * j 


> 


I \ -Bt  ! f 
uass  (c):  6 ■ = 


Xt 


L 0 


t€Xt 

_xt 


(2d). 


It  will  be  necessary  to  consider  several  cases  in  order  to  solve  the 
minimal  problem.  However,  there  are  some  results  that  can  be  established 
for  all  two-dimensional  systems.  In  this  connection,  let  T + be  the  curve 
obtained  by  starting  at  the  origin  and  following  the  P-curve  backwards  in 
time.  I*  + is  the  curve 


ds  J w 


(t  < 0) 


(29) 


Also,  let  F be  the  N-curve  symmetric  to  F+  with  respect  to  the  origin. 
The  following  theorem  will  now  be  proved: 

Theorem  3 ' In  the  two-dimensional  system 


x = Ax 


(e  f 0) 
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i 

i 
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where  A has  real  eigenvalues,  the  minimal  paths  for  points  on  T + T 
are  just  the  portions  of  F + or  T connecting  the  initial  point  and 
the  origin. 

Let  z°  be  any  point  on  T + s except  the  origin.  Let  the  minimal 
path  described  in  Theorem  3 be  denoted  by  P + and  let  S be  any  other 
path  from  z1"  „ Iu  must  be  shown  that  the  time  length  of  T 0 is  less 
than  the  time  length  of  5 . This  will  be  done  by  projecting  both  paths 
on  the  z^-axis.  It  will  turn  out  that  the  time  lengths  of  the  paths  are 
no  altered  by  projection.  It  will  further  turn  out  that  the  projected 
paths  are  paths  in  the  z^-axis  for  the  one-dimensional  minimal  problem 
which  was  considered  earlier.  From  (27),  the  equation  of  the  z^-component 
of  any  P-  or  N-curve  can  be  obtained,  the  result  being  the  same  regardless 
of  which  of  the  three  forms  the  matrix  D takes.  This  is  easily  seen  to 
be 


o + 2 \ 

« (Z  2 - T ) 


if  X / 0 and 


z 


2 


+ 


(30) 


(31 


\ 

/ 


I 

I 

t 

v 

r 

t- 

L 

l 

I 

I 

I 

I 


if  X - 0.  It  can  be  seen  that  the  time  length  of  a P-  or  an  N-arc  de- 
pends only  on  the  z,£  coordinates  and  is  not  altered  by  projection  on  the 
z^-axis.  It  is  also  seen  that  (30)  and  (31)  are  just  the  P-  and  N-curves 
of  the  one-dimensional  problem 


Z2 


(32) 

. i 


which  was  discussed  earlier.  For  the  one-dimensional  problem,  the  pro- 
jection of  r is  the  P-arc  connecting  z°„  with  the  origin,  which  is 
just  the  minimal  path  for  the  one-dimensional  problem.  The  projection  of 
S , which  must  contain  at  least  one  N-arc,  will  be  different  from  the 
unique  minimal  path  of  the  one-dimensional  problem.  The  time  length  of 

r +Q  is  therefore  less  than  the  time  length  of  S . This  completes  the 
z 

proof. 

[ 

. u 
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In  some  degenerate  cases,  the  minimal  paths  just  found  are  the  only 
ones.  These  cases  are  given  by  the  following  theorem: 

Theorem  1*:  If  e is  an  eigenvector  of  the  matrix  A , then  minimal 

paths  exist  only  for  points  on  T+  or  T . 

It  will  be  shown  that  the  P-curve  through  the  origin  is  a straight 
line  and  therefore  so  is  the  N-curve  through  the  origin  (by  symmetry). 
Since  the  N-curve  and  the  P-curve  coincide  and  since  the  N-  or  P-curves 
through  any  point  are  unique,  no  other  P-  or  N-curve  can  intersect  T 
or  T " . Thus,  there  is  no  way  to  reach  the  origin  from  a point  not  on 
F ’ or  F . 

In  terms  of  the  z- coordinates,  the  equation  of  the  P-curve  through 
the  origin  is 


z 


ds-w 


(33) 


Differentiating  (33)  gives 


z 


ds-Dw  ♦ w 


(32*) 


Clearly,  e is  an  eigenvector  of  A if  and  only  if  w is  an  eigenvector 
of  D . Therefore,  Dw  = Xw  and  (3U)  becomes 


z = X 6 


Dt 


/ £-Ds 
•'n 


ds-w  + w 


(35) 


or 

z = Xz  + w 

This  is  the  equation  of  a straight  line  and  the  proof  is  complete. 

The  eigenvectors  of  D , corresponding  to  ^ses  (a),  (b),  (c)  of 
(26),  are 
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rr 

j 

0 

Case  (a): 

1 0 . 

and 

1 

Case  (b): 

any 

vector 

) 

" 1 " 

Case  (c): 

0 

• 

1 

i 

i 

! 

r (36) 


According  to  the  theorem  just  proved,  case  (b)  is  a degenerate  case  ana 
need  net  be  considered  further.  To  avoid  degenerate  cases  in  (a)  and  (c), 
and  for  convenience,  it  is  assumed  in  the  sequel  that 


w.  > 0 
l 


(i  = 1,2) 


(37) 


Distinct  Nonzero  Eigenvalues 


In  this  case,  the  matrix  D assumes  a diagonal  form.  The  DE's 
are  in  the  uncoupled  form 

z = Dz  - w , (38) 

where 


rH 

CQ 

1 

i 

t-1 

o 

r— 

( — i 

1* 

L 

z = 

, D * 

> w = 

z 

0 X, 

w,. 

L 

2 

2 

X^  and  X^  are  the  eigenvalues  and  X^  / X^  . The  equations  of  the 
P-  and  N-curves  are 


+ Dt.  o + 

- u = £ ( z - 


u; 


(39) 


u = 


D~"V 


"wl" 

X1 

1 

J 

1 

€ " 0 

W2 

- Dt 

X t 

0 f 

r^l 

where 
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The  nature  of  the  P-  and  N-curves  depends  on  the  sign  of  x^  and 
Xg  • There  are  three  possibilities:  (l)  XX  0 , X^  < 0 , (2)  X.,  > 0 , Xg  > 0, 

and  (3)  X <C  0 , X^  > 0 - These  correspond  to  the  eases  where  the  singu- 
larity of  the  differential  equation  is  (l)  a stable  node,  (2)  an  unstable 
node,  and  (3)  a saddle  point.  Sketches  of  the  P-  and  N-curves  are  shown 
in  Figures  It,  5,  and  6,  pages  20  and  21.  These  Figures  are  drawn  under  the 
assumptions  that  w1  > 0 and  w^  > 0 . In  case  (l),  the  P-  and  N-curves 
are  semi-parabola-like  curves  that  tend  to  +u  as  t-*-*oo.  In  case  (2), 

the  curves  have  the  same  shape  but  with  reverse  direction.  The  reason  for 

+ 

the  interchange  of  -u  is  that  the  components  of  w are  assumed  to  be 
positive.  In  (3),  the  P-  and  N-curves  are  hyperbola-like  curves  around  the 
points  +u  . 

The  case  of  distinct  negative  eigenvalues  will  be  discussed  first. 

This  will  serve  as  a model  for  the  rest  of  the  proofs.  The  minimal  paths 

from  points  of  V and  T are  already  known.  It  is  necessary  to  find 

minimal  paths  from  all  the  other  points  in  the  plane. 

Let  C be  the  simple  curve  T ++  T - + 0 . This  curve  divides  the 
plane  into  two  open  regions.  The  region  which  contains  the  point  +u 


FIGuRS 
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will  be  called  K , and  the  region  which  contains  the  point  -u  wil'' 

be  called  R . If  z°  is  any  point  in  R , the  N-arc  from  z°  must 

intersect  T at  some  point  q , since  the  N-arc  connects  z°  and  +u 
which  are  on  opposite  sides  of  C . Also,  this  N-arc  cannot  intersect 
r since  r is  an  N-arc  and  the  N-arc  through  any  point  is  unique. 

Thus,  z°qO  is  a path  from  z°  (it  will  turn  out  that  this  is  the 

minimal  path).  A similar  argument  applied  to  points  of  R shows  that 
at  least  one  path  exists  from  every  point  in  the  plane.  The  points  in 
R can  also  be  obtained  by  starting  at  points  of  1 and  following 

the  N-arc  through  the  point  backwards  in  time,  that  is,  for  a time  t < 0. 

A similar  statement  holds  for  points  in  R*  . 

The  following  theorem  provides  the  minimal  paths: 

Theorem  5:  In  the  case  of  distinct  negative  eigenvalues,  a unique 

minimal  path  exists  from  every  point  in  the  plane.  If  z°  is  in  R , 
the  minimal  path  consists  of  followring  the  N-arc  through  z°  until  it 
intersects  T + and  then  following  T + to  the  origin.  The  minimal 
paths  for  points  in  R*  are  obtained  by  symmetry.  (In  terms  of  the  con- 
trol function  f , this  means  f is  -1  in  R and  on  T , and  f 
is  +1  in  R and  on  T . ) 

Let  z°  be  any  point  in  R . If  any  path  from  z°  leaves  R , 
it  must  pass  through  the  curve  T + . The  path  cannot  leave  R by  • 
passing  through  T since,  as  cai  easily  be  proved,  all  curves,  neces- 
sarily P-curves,  intersecting  P point  into  R . After  a path  in- 

’+  + 
tersects  I , the  minimal  path  from  that  point  must  coincide  with  1 

Therefore,  only  those  paths  from  z°  that  lie  entirely  in  R and  whose 
last  corner  is  on  P need  be  considered. 

Let  the  minimal  path  described  in  Theorem  5 be  called  A . It  is 
known  that  A exists  and  is  unique}  also,  A has  exactly  one  corner. 

A is  the  path  z°q0  in  Figure  7.  It  will  be  shown  that  the  time  length 
of  A is  less  than  any  path  that  contains  two  corners.  Successive  use 
of  this  result  will  prove  that  A has  shorter  length  than  any  other  path. 

Consider  any  path  from  z°  with  two  comers,  that  is,  for  a time 
To,  follow  a P-arc  from  z°  to  a point  z^  , then,  for  a time  , follow 
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ar  N-arc  to  a point  z^  (on  T )>  then  for  a time  T ^ , follow  T+  to 

the  origin.  If  the  point  z°  .ls  kept  fixed,  T can  be  taken  as  the  in- 

derendent  variable,  since  given  T , then  TL  and  T,  can  be  determined 

o i d. 

uniquely.  Let  the  total  time  length  of  this  path  be  I"  , that  is, 

r = rQ  *t1  +?2  T0’zi’T't  > o)  * (ij0) 

It  will  be  shown  that 


^>0  for  ro>0  . 


The  minimum  time  will  then  occur  when  T = 0 . However,  T = 0 gives 

o o 

the  desired  path  A . 

By  using  the  equations  of  the  P-  and  N-curves,  the  following  equa- 


i • 

tions  for 

o 1 

z , z , 

and  z^ 

■ - 

o 

-DT). 

z ' + 

c 

II 

A 

r\ 

u 

i > 

1 

-BT1 

(z  - u) 

i T 

z - 

u = € 

fi 

o 

z + 

-Dr 
, o 

u = € 

(z1  * u) 

By  successive  elimination  , it  is  found  that 

! -D(r  +r  +r  ) _d(t  +t  \ _dt  ] 
z + u = ^ c - c c + It  ju  (U3) 

or,  taking  components, 

f -\  v -v  . -7-  \ ■>  -r  i 

z°i  + Ui  = 1 6 1 ~2€  1 ° 1 +26  1 °|ui  (i  = 1,2)  . (hh) 


Differentia  ting  both  sides  with  respect  to  and  solving  for  d.T/dJ^  = T 

gives,  after  some  simplification. 


• xir2  • MW 

T-  2€  1 d (r  + 1)  = -26  1 1 * 


(i  = 1,2)  . (U5) 
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These  are  two  simultaneous  equations  for  T and  (T  + l)  . Solving  for 
r by  determinants  results  in 


( VVt2 

i ^ ^ 


X,T 


2 € 


1 1 


1 € 


X2*L1 


1 € 


i e 


X1T2 


V* 


(U6) 


The  determinant  in  the  denominator  is  not  zero  since  X^  / X^  and 

> 0 . Clearly,  the  quotient  of  the  determinants  will  be  positive  since 
the  numerator  will  have  the  same  sign  as  the  denominator.  It  has  thus 
been  shown  that  t > 0 for  Tq  > 0 , and  the  proof  is  complete. 

The  case  of  positive  distinct  eigenvalues  can  now  be  disposed  of 
fairly  readily.  A study  of  Figure  5 suggests  that  no  paths  exist  from 
points  outside  the  open  region  Q bounded  by  the  P-curve  through  +u 
(followed  backwards)  and  the  N-curve  through  -u  (followed  backwards)  as 
shown  below: 

Z2 


VTCiinv  A 

Of  course  this  can  be  proved  analytically  (see  Reference  l)  and  is  based 
on  the  fact  that  all  curves  (necessarily  N~curves)  intersecting  the  P-arc 
connecting  -u  and  +u  are  directed  cut  of  the  region  Q . 
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The  proof  given  for  Theorem  5 holds  good  without  change  for  points 
in  Q since  the  only  place  that  the  sign  of  the  eigenvalues  came  into 
the  proof  was  in  guaranteeing  the  existence  of  paths.  The  region  R is 
now  the  region  in  Q and  between  the  curve  C=  r+  * T +0  and  the 
N-curve  through  -u  , and  R is  the  reflection  of  R in  the  origin. 
The  following  theorem  therefore  holds  (Figure  d shows  a typical  minimal 
path) : 

Theorem  6:  In  the  case  of  positive  distinct  eigenvalues,  a minimal 

path  exists  (and  is  unique)  only  in  the  region  Q defined  above.  If 
the  regions  R+  and  R are  redefined  as  above,  the  minimal  paths  are 
the  same  as  in  Theorem  5- 

In  the  case  where  X,  < 0 and  X.p  > 0 , the  only  changes  necessary 
I'cr  the  above  theorem  to  remain  true  are  in  the  definitions  of  Q , R , 

R . Q is  easily  determined  to  be  all  points  in  the  horizontal  sv.-;p 
—Up  < Zp  < Up  . The  region  R is  the  region  between  the  curve  C and 

the  line  = Up  , and  R , its  reflection  in  the  origin.  A typical 
minimal  path  is  shown  in  Figure  9- 


-u 


FIGURE  9 


Distinct  Eigenvalues  with  One  Eigenvalue  Zero 

The  matrix  i>  still  assumes  the  diagonal  form 


0 

0 


4 0) 
* / ' » 


<\>n\ 
\ *+  • / 


tsnrmwm  *wf$wtvw^7^w  tywarw  ?«»t^ 


-26- 


The  solutions  of  the  DE  z = Dz  - w are  different  from  those  of  the  last 
section  since  d is  now  a singular  matrix.  The  solutions  are 


+ W1  _ -Xt I o ♦ wl\ 

zl  " x c \2  l ~ TJ 


o + . 

z2'  z Z - V ■ 


The  sketches  of  these  P-  and  N-curves  are  shown  in  Figures  10  and  11, 
page  27,  for  the  cases  X < 0 and  X > 0 , respectively. 

It  can  be  seen  that  if  X < 0 , a path  exists  from  every  point  in 
the  plane  while  if  X>  0 , paths  exist  only  from  the  points  | z-^  j < |wj/x| 
The  same  method  of  proof  will  be  used  as  in  the  last  section.  The  theorem 
to  be  proved  will  not  be  reformulated  since  it  is  the  same  as  Theorem  6, 
where  Q is  the  region  from  which  paths  exist  as  defined  above,  and  R 
and  R+  are  the  regions  above  and  below  the  curve  C , respectively. 

The  notation  of  Theorem  5 will  be  used.  It  is  only  necessary  to 
find  the  equation  for  T based  on  the  new  P-  and  N-curves  and  see  if 
T>  0.  Since  the  first  equation  of  (U6)  is  the  same  as  the  z.. -component 
of  (39) j one  equation  for  X is  ()*5): 


T-  26  * (T  + 1)  - -2€ 


The  second  equation  can  easily  be  found  by  using  the  second  equation  of 
(ho).  Thb  equations  similar  to  (1*2)  are 


z 2 = v2^~X2^ 


z12  = z22  - ^(-V 


*°2  * l\  * v2(-To)  ' 


By  elimination,  it  is  found  tnat 


i^w**««*  «*•*«*?**<  W»*W  g*”^1  **»■•; 
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O — — t -T* 
2 2 ' "2V  “2 


which  is  differentiated  with  respect  to  T to  obtain 

o 


t1  - r2  - l = o 


It  is  recalled  that  T = T + + T 2.  * differentiating  this  equation  and 

combining  it  with  (52)  gives 


.1  r 

"l  2 


Substituting  this  into  (u9)  and  solving  for  T results  in 


r=  26 


VC 

.2 1-6  1 


1-6 


Sirce  X / 0 and  > 0 , it  is  seen  that  ^ > 0 for  ro  > 0 , 

which  completes  the  proof. 


Repeated  Eigenvalues 


The  matrix  D now  takes  the  form 


r>.  r 

D = ! 

L°  x_ 

If  X / 0 , the  P-  and  N-curves  are 


+ Dt  . o + \ 

z - u = 6 (z  - u)  , 


where 


r,-l 

u = D w 


Dt  = \t 


or  if  X = 0 , the  equations  are 


R-h$9 
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+ t . . / o + s . o 

’•1  “2  T tv*  * * *XJ  ’ * 1 


+ . . 0 
Z2  = “w2b  z 2 


;V  (5v) 


The  sketches  of  these  corves  are  shown  in  Figures  12,  13,  1)4,  pages 
30  and  31  o It  can  be  seen  that,  if  X <.0  , paths  exist  from  every  point 
in  the  plane,  while  if  X > 0 , paths  exist  only  in  a region  Q bounded 
by  the  P-curve  through  +u  followed  backwards  and  the  N-curve  through 
-u  followed  backwards. 

The  region  R is  the  region  above  the  curve  C if  X <0  and 
the  region  below  C if  X > 0 . With  this  notation,  Theorem  6 remains 
true.  The  proof  is  essentially  the  same  as  that  of  Theorem  6 and  will  not 
be  given.  The  only  change  is  in  the  algebra  due  to  the  different  equations 
for  the  P-  and  N -curves. 


Summary  of  Two-Dimensional  Cases 

Some  definitions  are  reviewed  below  for  convenience  in  summarizing 
the  results  of  the  two-dimensional  cases: 

Let  r + be  the  point  set  obtained  by  following  the  P-curve  through 
the  origin  for  all  time  t < 0 . 

Let  R be  the  point  set  obtained  by  starting  at  each  point  of  T + 
and  following  the  N-curve  through  the  point  for  all  time  t < 0 . 

Let  r and  R+  be  the  point  sets  obtained  by  reflecting  H + 
and  R in  the  origin. 

Let  C = r + + T“  + 0 . 

Let  Q = R + R + C . 

The  following  theorem  summarizes  the  results: 

Theorem  7:  In  the  two-dimensional  sysrem 

x = Ax  +ef(x)  , 


where  e / 0 , f = -1  , and  A has  real  eigenvalues,  a unique  minimal 


R-U59 

-32- 


path  exists  for  points  in  Q , and  no  paths  exist  for  points  outside  of 
Q . Jot  points  in  R , the  minimal  path  consists  of  following  the  N-arc 
through  the  point  until  it  intersects  P+  and  then  following  T + to 
the  origin.  Minimal  paths  for  points  in  R+  are  obtained  by  symmetry. 

(In  terms  of  the  control  function  f , this  means  that  f = -1  in 
R + r and  f = +1  in  R++  T+.) 

The  curve  C is  the  common  boundary  of  the  regions  R+  and  R~  .. 

It  is  on  this  curve  alone  that  the  minimizing  control  function  changes 
sign  or  "switches"  from  the  value  +1  to  the  value  -1  or  vice  versa. 

For  this  reason,  C is  called  the  switching  curve. 

The  region  Q from  which  minimal  paths  exist  depends  on  the  nature 
of  the  eigenvalues.  In  the  degenerate  case  where  the  vector  e is  an 
eigenvector  of  A , the  region  Q degenerates  into  the  curve  C . Aside 
from  this  case,  Q is  a "two-dimensional"  region.  If  the  eigenvalues 
are  nonpositive,  Q is  the  entire  plane.  If  both  eigenvalues  are  posi- 
tive,  Q is  a bounded  portion  of  the  plane,  the  boundary  curves  being 
the  P-curve  through  +A  1e  and  the  N-curve  through  -A-1e  (both  followed 
backwards).  If  only  one  eigenvalue  is  positive  while  the  other  is  non- 
positive, Q is  an  unbounded  portion  of  the  plane  between  two 
lines,  these  straight  lines  being  solutions  of  the  DE  that  are  parallel 
to  the  eigenvector  corresponding  to  the  positive  eigenvalue. 

It  should  be  remarked  here  that,  in  an  (uncontrolled)  physical 
system,  positive  eigenvalues  mean  some  sort  of  instability.  However,  it 
has  been  shown  above  that,  by  the  proper  use  of  an  "on-off"  control,  the 
system  can  be  made  stable,  provided  the  error  and  its  derivative  remain 
in  the  region  Q . 

Also,  another  way  of  characterizing  the  minimal  paths  is  as  follows: 

If  any  paths  at  all  exist  from  a point,  then  a unique  minimal  path 
exists.  The  minimal  path  is  the  path  with  a minimum  number  of  corners. 

The  Minimal  Problem  for  x^  = f(x^,x^)  with  j f | < 1 

The  broader  minimal  problem  will  be  considered  for  the  simple  DE 

x^  = f(x1,x1) 

* Includes  repeated  eigenvalues. 
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or  the  equivalent  system 


Xr,  ~ f(X-,,X„) 

C.  X L. 

where  | f(x. ,x^)  j < 1 . It  will  be  shown  that  the  minimizing  function  is 
still  the  sane  function  as  the  one  found  among  the  restricted  class  of 
functions  capable  of  taking  on  only  the  values  *1  . 

In  order  to  guarantee  solutions  to  (5b),  it  is  necessary  to  put  some 
restrictions  on  f ; however,  f is  still  allowed  to  have  some  discon- 
tinuities so  that  it  will  include  admissible  functions  of  the  type  f = -1  . 
It  is  recalled  that  in  the  case  f = -1  , the  only  admissible  functions 
were  those  such  that  the  solution  curve  from  any  point  consisted  of  a finite 
number  of  alternating  P-  and  N-arcs.  This  is  equivalent  to  the  requirement 
that  any  bounded  region  can  be  split  up  into  a finite  number  of  subregions, 
in  each  of  which  f is  continuous. 

In  the  broader  problem  (58),  it  is  assumed  that  any  bounded  region 
can  be  split  up  into  a finite  number  of  subregions,  in  each  of  which  f is 
continuous  and  satisfies  a Linshitz  condition.  However,  f should  be  de- 
fined at  every  point  in  the  plane.  This  assures  a unique  solution  from 
every  point  in  the  plana  if  it  is  required  that  the  solution  curve  be  a 
continuous  curve  (in  the  phase  plane). 

Let  A be  the  minimal  path  obtained  previously  when  f = -1  , and 
let  ji  be  any  solution  to  (5o)  with  the  same  initial  point  as  A . Penote 
T as  the  time  length  of  a solution  curve.  With  this  notation,  it  is  only 
necessary  to  prove  the  following  theorems 

Theorem  0;  I"  (A)  S~(h)  . 

It  will  be  sufficient  to  prove  the  theorem  for  initial  points  in 
R~  and  above  the  x, -axis.  The  path  A from  such  an  initial  point  x° 

X 

is  shown  on  the  following  page  (the  path  x qrO): 


) 

> (50) 
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An  inequality  for  any  arc  of  n will  now  be  obtained.  Since 

-1  < x2  < +1  . (59) 

Integrating  this  from  0 to  t gives 

-t  < Xg  “ r-J 2 — **  (60) 


M > 


(61) 


where  the  equality  sign  holds  only  for  a P-  or  an  N-arc.  I'his  expression 
proves  the  theorem  immediately  for  points  oh  T+  or  T »- 


An  inequality  for  the  slope  of  n is  obtained  next;  Let  the  slope 
of  a be  m ; then 


VUC./ 


I 


i 

\ 

I 

! 

i 

* 

i 


or 


j 


f 


I 


(63) 


Therefore,  the  slope  of  n at  any  point  lies  between  the  slope  of  the 
N-curve  and  the  si ore  of  the  P-curve  through  the  point. 

Consider  the  curve  n through  the  point  x°  „ Because  of  (63),  n 
cannot  pass  through  either  the  P-curve  through  x°  or  the  N-arc  x°q 
through  x°  . If  n is  to  reach  the  origin,  it  must  therefore  cut  the 
x,-axis  at  a point  q1  to  the  rigrt  of  (or  possibly  at)  the  point  q. 

i. 

From  (6l),  it  is  seen  that  T (x°q)  < I"(x°q  ) . 

Now  consider  the  rest  of  n . If  n passes  below  the  horizontal 
line  through  r , (61)  shows  that  T (A)  < T(ji)  , and  the  theorem  is 

proved.  Otherwise,  if  n is  to  reach  the  origin,  it  must  intersect  A 
at  some  point  s along  the  arc  qr  . For  this  portion  of  n , it  is 
seen  again  from  (6l)  that  T (x°qs)  <,r(x°q^s)  . It  is  only  necessary  to 
prove  a similar  inequality  for  the  remaining  part  of  n . There  are 
three  possibilities  to  be  considered: 

1.  n reaches  the  origin  by  staying  entirely  within  the  area 
bounded  by  the  part  of  A}  qrO  s and  the  x-^-axis. 

2.  it  leaves  this  region  through  the  arc  rO  . 

3.  ti  leaves  this  region  through  the  x^-axis. 

In  case  (l),  T (A)  < r(n)  since  the  time  length  along  any  solution 
is 


and  j x^  j along  n is  not  greater  than  j x^j  along  A. 

In  case  (2),  suppose  that  n leaves  the  region  for  the  first  time 
at  the  point  v on  the  arc  rO  . Then  the  desired  inequality  will  hold 
up  to  v by  case  (l)  and  after  v by  equation  (61)  so  that  T (A)<  r(n) 
will  still  hold. 
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I 


I 

\ 

I 


In  case  (3),  as  soon  as  r,  gets  to  a point  above  the  x^-axis,  the 
considerations  discussed  earlier  will  prevail  and  n will  have  to  in- 
tersect the  x^-axis  again  to  the  right  oi’  its  first  intersection.  Thus, 
the  length  of  n has  been  made  longer  before  case  (l)  or  (2)  applies, 
and  the  desired  result  follows. 

It  is  fairly  clear  that  the-  strict  inequality  will  hold  unless  A 
is  identical  with  n , but  this  means  that  f would  be  just  the  mini- 
mizing function  previously  found.  This  completes  the  proof. 


I 

I 

■ 

1 

1 

I 
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N -DIMENSIONAL  SYSTEMS 


Introduction 


Consider  the  minimal  problem  for  the  n-dimensional  system 


x = Ax  - e 


(e  i 0) 


(6k) 


where  A is  an  n by  n matrix  with  real  eigenvalues.  By  an  appro- 
priate nonsingular  substitution  x = Tz  , (6ii)  can  be  reduced  to 


z = nz  - w 


(w  / 0) 


(65) 


where  w = T ^e  and  D = T ^AT  , and  D is  in  Jordan  canonical  form, 
that  is. 


D = 


A-, 


a2  U 


O 


(66) 


Each  of  the  Ai  i3  a matrix  having  one  of  the  following  forms: 


X,. 

c 


o 


Case  (a): 


o 


o 


Case  (b): 


o 


(67) 


H~u5^ 

-.38- 


X 1 
X 1 


Casa  (c): 


X 1 
X 


Here,  X^  and  X are  eigenvalues  of  A . 

The  equations  of  the  P-  and  N-curves  are  still  given  by 

*t 

Dt  o + Dt  / -Ds  , 

Z = £ z - £ J £ ds-w 

*o 

The  matrix  function  6 ^ is 


O 

O''-. 


A ^ 

where  € has  one  of  the  three  forms  corresponding  to  cases  (a),  (b), 
(c)  of  (67): 


* o 


Case  (a): 


u 


x.t 


Case  (b)t  g i , 
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1 t 


21 


x t 


t”’1 
(n-l) ' 

tn  2 


,,  , x xt ! 

case  (o):  £ 


1 t 


V (70) 


The  minimal  problem  will  be  solved  first  for  the  cases  where  D 
is  equal  to  one  of  the  matrices  in  (6?),  and  then  for  the  general  case. 
Certain  exceptional  cases  can  occur  ard  will  be  disposed  of  in  the  next 
section. 


Degenerate  Cases 

As  was  seen  in  the  two-dimensional  case,  certain  degenerate  cases 
can  occur.  The  following  theorem  carries  over  directly  from  the  two- 
dimensional  case. 

Theorem  9'  If  w is  an  eigenvector  of  the  matrix  D , then 
minimal  paths  can  exist  only  for  points  on  T + + f . 

I and  r are  defined  just  as  in  the  two-dimensional  case.  The 
former  proof  holds  good  without  change.  It  is  recalled  that  T+  and  T 
lie  on  a common  straight  line  through  the  origin. 

Even  if  w is  not  an  eigenvector  of  D , there  still  can  occur 
degenerate  cases.  Suppose,  for  instance,  that,  in  the  system  z * Dz  - w , 


^1  1 0 ^ 

2 

H O 

1 

0 

D = 

o 

and  w “ 

a 

o 

.O  ’”Vx_ 

o 

0 

s 


(71) 
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where  the  X .....  X , are  distinct.  It  is  easily  seen  that  w is 
1 n-1 

not  an  eigenvector  of  D . However,  consider  the  system 

- =_  + _ 
z = Dz  - w , 

where  z denotes  the  n-1  dimensional  vector  whose  components  are 

,z  i and 
c 5 n 


It  can  be  seen  that  w is  an  eigenvector  of  D corresponding  to  the 
eigenvalue  . In  this  n-1  dimensional  system,  the  equations  of 
the  P-  and  N- curves  are  obtained  simply  by  leaving  off  the  equation  for 
z.^  in  the  n- dimensional  system.  If  a path  from  z in  the  n-1  di- 
mensional system  is  considered,  a unique  path  to  a unique  point  z in 
the  n-dimensional  system  can  be  determined  simply  by  working  backwards 
frem  the  origin  and  following  P-  and  N-curves  for  the  same  lengths  of 
time  as  for  the  path  in  the  n-1  dimens icral  system.  Conversely,  every 
path  from  z will,  in  the  same  manner,  determine  a unique  path  from  its 
projection  z . If  no  path  exists  from  z , no  path  will  exist  from  z . 


Now,  according  to  Theorem  9,  paths  can  exist  in  the  n-1  dimensional 
system  only  for  points  on  a certain  straight  line  through  the  origin. 
Therefore,  paths  can  exist  in  the  n-dimensional  space  only  for  points 
which  project  onto  this  straight  line.  These  are  points  with  z ^ arbi- 
trary and  Zg  > , •••  , zn  satisfying  an  equation  of  the  form 


Z2 


W2 


w 


2 


Thus,  minimal  paths  can  exist  only  for  a two-dimensional  plane  through 
the  origin. 


l 


I 

I 

! 

i 

I 

I 

I 

I 


T 

l 

' f 


5 


I 
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lt is  pointed  out  that  the  above  demonstration  depended  on  the  fact 
that  the  equations  for  * "*  > zn  did  not  contain  z^  . If  z^ 

were  contained,  then  z would  not  be  a solution  of  the  same  type  of  dif- 
ferential equation  as  z . 

Reasoning  si  mi jar  to  the  above  will  prove: 

Theorem  10:  Suppose  the  k-dimensional  system 

z = Dz  - w (W  j/  0) 

can  be  formed  bv  taking  k of  the  component  equations  of  the  n-dimen- 
sicnal  syst  em 

z = Dz  - w (w  / 0) 

If  w is  an  eigenvector  of  D , then  minimal  paths  can  exist  only  in  a_ 
certain  n-k+1  dimensional  hyperplane  through  the  origin. 

In  order  to  avoid  these  exceptional  cases,  it  is  assumed  in  the 
sequel  that  (l)  if  X is  a repeated  eigenvalue,  X appears  only  in  a 
diagonal  matrix  of  the  form  (c)  of  (67)  and  in  only  one  such  matrix, 
and  (2)  w\  / 0 for  i = 1,2,  , n . Such  a system  is  called  a 

nondegenerate  system.  Theorem  10  holds  only  for  k*i  for  a nonde- 
generate system,  and,  in  fact,  it  will  be  seen  that  minimal  paths  will 
exist  for  all  points  in  some  n-dimensional  neighborhood  of  the  origin. 

Terminology  and  Statement  of  the  Main  Theorem 

In  finding  the  minimal  paths,  it  will  again  be  convenient  to 
work  backwards  from  the  origin.  Using  this  procedure,  the  following 
point  sets  or  "surfaces"  in  n-dimensional  space  are  defined: 

IL.  s The  point  set  obtained  by  following  the  P-curve  through  the 
origin  backwards  in  time.  R^  is  the  set  of  points  z satisfying 

-ti 

z^  = € ■LJ*  € ds-w  (t.  >0) 

o 1 
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R^s  The  point  set  obtained  by  starting  at  any  p CiTit  of  and 

following  the  N-curve  through  the  point  backwards  in  time,  that  is,  points 
2 

z such  that 


2 

z 


ds-w 


(tk>  0)  . 


Rn:  The  points  z satisfying 


-Dt  , , 

n ^ n n-1  . *n-l 

z = c z + (-1)  € 


-Dt  ftr 

■ "I 


-Ds  . 

6 ds.w 


(t  >0)  . 


-1 


Let  R^  be  the  point  set  obtained  by  reflecting  R^  in  the 


origin. 


ut  Qk  ’ \ * \ 


-1 


Let  Q*  - 


i • 


i*l 

R^  is  the  same  as  the  curve  T + previously  defined  and  R^  1 

is  the  same  as  T ~ . consists  of  all  points  from  which  a path  exists 

with  no  corners;  consists  of  all  points  from  which  a path  exists 

with  exactly  cue  corner;  and,  in  general,  Q,  consists  of  all  points  from 

* «#• 

which  a path  exists  with  exactly  k-1  corners.  is  the  logical  sum 

of  the  sets  and  consists  of  all  points  from  which  paths 

exist  with  k-1  or  less  corners. 


Two  theorems  will  now  be  formulated.  To  prove  these  theorems,  it 
is  necessary  to  consider  several  cases  which  will  be  done  in  the  following 
sections. 

Theorem  11:  In  a nondegenerate  system  with  real  eigenvalues,  the 

point  sets  y • • - ,Qn  are  mutually  exclusive.  From  each  point  in 

, a unique  path  exists  with  exactly  k-1  corners,  and  no  path  exists 
with  less  than  k-1  corners. 


i 
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The  main  theorem  of  this  report  is: 

Theorem  12:  In  a nondegenerate  system  with  real  eigenvalues,  a_ 

unique  minimal  path  exists  from  points  in_  the  region  Qn  , and  no  paths 
exist  outside  of  . The  minimal  path  from  a point  iri  is^  the 

unique  path  with  exactly  k-1  corners. 

As  is  seen  from  the  above  theorems,  all  minimal  paths  contain 
n-1  or  less  corners.  All  the  corners  of  minimal  paths  lie  on  the  sur- 

-X- 

face  Q ^ , which  is  called  the  switching  surface. 

Distinct  Eigenvalues 

Consider  the  case  where  the  matrix  D in  the  nE 
z = Dz  - w (wm  / 0)  (73) 

has  the  form 


where  the  eigenvalues  X^  , ...  , X^  are  distinct.  It  will  be  assumed, 
for  the  present,  that  no  eigenvalue  is  zero.  The  equations  of  the  P-  and 
N-curves  are  then 


+ Dt  , o + x 
z - u =€  (z  - u) 


where  u = D and  u.  / 0 since  w.  ^ 0 . 


j.  i 

Theorem  13:  The  sets  , . . . , are  mutually 

From  each  point  in  , a unique  path  exists  with  exactly 
and  no  path  exists  with  less  than  k-1  corners.  No  paths 

-H- 

points  outside  of  . 

This  theorem  will  be  proved  by  induction.  The  proof 
given  earlier.  The  idea  of  projection  will  be  used  again. 


(75) 


exclusive. 
k-1  corners, 
exist  from 


for  n •-  2 was 
Denote  the 


n-1  dimensional  system  obtained  by  leaving  off  the  first  equation  in  (73) 

by 


z = Dz  ± w . (76) 

^ denotes  the  regions  in  the  n-1  dimensional  space  with  which  the 
theorem  is  concerned.  By  the  inductive  assumption,  the  theorem  holds 
for  , . . . , ^ . - Consider  any  point  z in  . A unique  path 

with  exactly  k-1  corners  exists  from  this  point.  This  path  consists 
of  a definite  sequence  of  P-  and  N-curves  followed  for  definite  times  in 
the  n-1  dimensional  space.  By  starting  at  the  origin  and  following  the 
same  sequence  backwards  for  the  same  times,  but  now  in  n dimensions,  a 
unique  point  z is  reached  by  means  of  a path  with  exactly  k-1  corners. 
This  path  is  unique  since  otherwise  the  corresponding  path  from  z would 
not  be  unique.  There  is  no  path  from  z with  less  than  k-1  corners 
since  no  path  from  z contains  less  than  k-1  corners,  z is  therefore 
in  . Similar  reasoning  shows  that  to  each  point  z in  there  is 

a unique  point  z in  ; therefore,  a unique  path  of  exactly  k-1 
corners  exists  frcm  z and  no  path  with  less  than  k-1  corners.  This 
proves  the  theorem  for  the  sets  Q.,  , , ...  , Qn  ^ . Incidentally,  it 

can  be  seen  from  the  above  that  a one-to-one  correspondence  exists  between 
Of  , Q2  , ...  , Qn_1  and  Q1  , , ...  , Qn_1  . The  switching  surface 

• 

Ci  • * 0 

i-1 

can  therefore  be  represented  as  a single-valued  function 

Z1  = F(Z2  > z3  > *..  » * 

From  the  definition  of  the  Q-sets,  it  also  follows  that  the  switching 
surface  is  continuous. 

The  theorem  for  the  set  Qn  has  not  yet  been  proved.  To  do  this, 

it  is  necessary  to  have  the  explicit  expression  for  points  in  Q . It  is 

-1  n 
recalled  that  Q = R + R . The  equation  for  points  in  R is  ob- 
n n n n n 
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tained  frcm  the  definition  of  the  R-sets  and  by  the  process  of  successive 
elimination.  If  zn  represents  any  point  in  Rn  , then 


_ _ j -f'(T-,  + T'+.  • -+T_) 

n . , -.n-i  ' _ L / n 

z + (-1J  u = l 6 


-D(r.  +...  + r)  , -DT  I 

-2€  * n + ...  + (-l)n_12fc*  nj  a , (77) 


where  all  the  T ^ are  greater  than  zero.  The  equation  for  points  in 
^ is  obtained  by  replacing  zn  by  -zn  . rarameiers  are  changed  by 
the  nonsingular  substitution 


n-l 


'N 


* T1  + x 2 * ♦ • • + T 

= TV, ' + . . • + T , 

C j 


n 


T , + T 
n-l  n 


(7d) 


J 


i ; 


It  is  noted  that  a^  > a^  > ...  > an  >0.  Equation  (77)  now  becomes 

zn  + (-1)"“^  = ' 6 1 - 2€  2 + ...  + (-l)”’3^  n ) u . (79) 

v.  j 

Now  if  zn  is  a point  in  Qn_-j_  > equations  (79)  are  satisfied  by 

certain  values  of  a.  . a,,  ...  ,a  . and  a = 0 . Consider  the 

1 ’ 2 n-l  n 

Jacobian  of  this  system: 


dJi 

aa 


j i 


k 


\„a 
€ 11 

^Xla2  c Xl3n 

€ X2ai 

\2&2  \2&n 

C # • c C 

0 

0 0 

• 

• 0 

0 

• • 

X a. 

s n 1 

X a0  X a 

n 2 _ n n 

C » 9 • C 

(SO) 


7 * 
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nore, K is  a nonzero  constant.  The  X.  are  all  different  from  aero  by 

x 

assumption,  and  the  u^  arc  all  different  from  zero  since  the  system  is 

nondegenerate.  The  determinant  is  also  different  from  zero  if  the 

and  the  a.  are  distinct,  which  is  the  case  here.  (The  determinant  is 
1 

a generalization  of  the  Vandermonde  determinant,  and  a proof  that  it  does 
not  vanish  is  given  in  Reference  3,  part  7,  problem  26.)  Equation  (7?) 
can  therefore  be  solved  uniquely  for  ai  - a2  * •••  , an  scme  neighbor- 
hood of  zn  . This  shows  that,  in  some  neighborhood  of  z ‘ , a unique 
path  with  exactly  n-1  comers  exists.  As  a is  increased,  the 
Jacobian  remains  different  from  zero  so  that  a unique  path  with  exactly 
n-1  corners  exists  from  every  point  in  . With  the  same  type  cf 
reasoning,  it  can  oe  established  that  no  path  exists  from  a point  in 
with  less  than  n-1  comers. 

To  complete  the  proof,  it  must  be  shown  that  no  paths  exist  from 
points  outside  of  . Consider  such  a point  z . If  a path  exists 
from  z , it  must  contain  more  than  n-1  corners;  say  it  contains  m-1 
corners  with  m > n . An  equation  similar  to  (79)  could  be  set  up  con- 
taining m parameters.  These  equations  could  always  be  solved  for  the 
first  n of  these  parameters  and  the  other  parameters  set  equal  to  zero. 

Therefore,  a path  exists  from  z with  n-1  or  less  corners  and  z must 

* 

therefore  be  in  . This  completes  the  proof. 

* 

A word  should  be  said  here  about  the  extent  of  , the  region 

of  existence  of  paths.  From  (79),  it  is  seen  that,  if  some  of  the  eigen- 

* 

values  are  negative,  the  region  Q will  extend  to  infinity  in  some  di- 

n 

rections,  while  if  the  eigenvalues  are  all  positive,  the  region  will 
be  bounded. 

More  precise  information  is  possible.  It  could  be  proved  that,  if 
the  eigenvalues  are  all  negative,  paths  exist  from  all  points,  while  if 
the  eigenvalue  is  positivo,  paths  do  not  exist  for  points  | z-t  | £.  | WjA-J  • 

In  the  latter  case,  paths  do  not  necessarily  exist  from  all  points  satis- 
fying | z±  | < | WjAJ  • 

" If  it  happened  that  this  solution  gave  a =0  for  a point  in  the  neigh- 
borhood of  zn  , then  the  point  would  certainly  be  in  Q ^ and  could 
not  be  in  because  of  the  uniqueness  of  the  solution. 


I 

t 
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I 

I 

I 
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I 

I 
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■ 

I 
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The  main  theorem  for  this  case  will  now  be  proved. 

Theorem  lit?  if  the  eigenvalues  are  distinct  and  not  zero,  a unique 
minimal  path  exists  for  all  points  in  . The  minimal  path  for  a point 
in  is  the  unique  path  that  contains  k-1  corners. 

This  theorem  will  be  proved  by  induction.  The  theorem  has  been 
proved  for  n = 2 „ 

The  theorem  will  be  proved  first  for  points  on  the  switching  sur- 
face Q*  . Consider  any  pcint  z in  for  k < n-1  and  its  pro- 
jection z in  as  in  the  proof  of  the  preceding  theorem.  The  unique 

path  A ^ ^ with  k-1  corners  from  z projects  onto  the  unique  path 
A^  ^ with  k-1  corners  fresn  z . The  paths  A and  A have  the 
same  time  length.*  In  the  n-1  dimensional  space,  A ^ ^ is  the  minimal 
path  by  assumption.  Any  path  from  z other  than  A would  project 
into  a path  different  from  A and  would  therefore  have  a longer  time 
length  than  A^  . Therefore,  A is  the  minimal  path  and  the  theorem 
is  proved  for  all  points  in  • 

It  remains  to  prove  the  theorem  for  points  in  Qn  . Suppose  z is 

anv  point  in  Q . It  is  known  that  there  exists  a unique  path  A , 
n n-1 

from  z containing  exactly  n-1  corners;  it  is  desired  to  show  that 

A . is  the  minimal  path. 
n-1 

Consider  any  path  at  all  from  z . This  path  will  have  a certain 

number  of  corners  and  will  then  intersect  Q , . Fran  this  point,  the 

n-1 

minimal  path  is  uniquely  determined.  It  is  only  necessary  to  determine 

the  behavior  of  the  minimal  path  fran  z until  its  first  intersection 

n^l 

with  Q , . For  convenience,  label  the  ooint  z as  z . Consider 
n-1  n+1 

the  path  A from  z that  has  exactly  one  corner  before  reaching 

Q ..  . In  other  words,  A is  the  path  obtained  by  starting  at  z"  * 
n-l  n jj 

and  following,  say,  a P-arc  for  a time  T,  > 0 to  a point  z , not  in 

il  ^ 

Q . , and  then  an  N-arc,  for  a time  T , S 0 to  a point  z , which 
n-1  n-1 

is  in  ^ . From  this  point,  the  minimal  path  is  determined.  It  will 

be  shown  that  the  time  length  of  A can  be  reduced  by  making  T smaller, 

n n 

* Any  path  and  its  "projected"  path  have  the  same  time  length  by  construction. 
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and  that  it  is  a minimum  if  T^  = 0 , This  will  show  that  the  time  length 

of  A , is  less  than  tie  time  length  of  A • Successive  use  of  this 
n-l  n 

result  will  snow  that  A , is  shorter  than  a oath  with  any  number  of  cor- 
n-1 

ners . 

Using  the  equations  of  the  P-  and.  N-curves,  the  equations  for  the 

successive  corners  of  A are  found  to  be 

n 


n+1  , . 

z + u =6 


-»(*)  n 

n (zn  + u) 


n 

z - u 


-Dr  , 

_ _ n-l/  n-l  v 
= € (z  - u) 


n-l 
z ■» 


-nr 

. n-2/  n-2 

= € (z  + u) 


1 . / , \n+l  DV  o . / , \n*l  , 

z + '-1)  u = c '^z  + (-1)  u . 

By  successive  elimination,  it  is  found  that 


(z°  = 0,  T±>  0) 

J 


f -DT  -D(T +T  ) . -D(T+...+T) 

t"' 1+  c-x)n' fl«  -[2f  n-2«  n n-1  ♦ ...♦  !-i)n_1te  1 

„ -D(r  +...+? ) 

+ (-l)n€  ° n u . 


Now  let 


To  * T1  * 


. . ♦ T 


T , + t 

1 n-l  n 


where  a&  represent?  the  total  time  length  of  A 
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Substituting  (63)  into  (d2)  gives 


n+1  , . ,n+l  ( /Dan  "DVi.  , _,n-l  *Dal 

z + v-x;  u = <;fc  -26  + ...+  (-1)  2c 


-Da 

M-l)n€  °fu 


(8U) 


or,  taking  components,, 


-X.a  1 


(.1)ne'Ai“=jUi  . m 


Let  T *»  an  be  the  independent  variable,  then  differentiating 
both  sides  of  (8$)  with  respect  to  results  in 

—X. . a — X . a , . — X . a 

0 = 2£  1 n - 2€  1 n"x  a + ...  + (-l)n_12<~  1 1 a, 

n-1  i 


_ -X.a 

+ (-l)fc  1 0 


(i  = l,...,n)  , (86) 


i | 


i 1 ? 

I i 


i { 

I 5 

! J 


where  a.  = da. /da  . These  are  n simultaneous  equations  for 
x x n 


ao  ’*1*  •"  ’Vl 


-X  a 
. 1 n 


Solving  by  determinants  for  a gives 

o 


a 

o 

T 


-X.a  . 
. 1 n-1 


-X  a 
- n ri 


-X  a 

- n n-1 


-X.a  . 
s 1 n-1 


-X  a . 
i ^ n n-1 


"X1S1 


o o 


. € 


-X  a. 
n 1 


-X.a 
. 1 o 


• • 


-X  a 
c no 

• • " 


(87) 


i i 


R -h59 
- 50  - 


The  determinant  in  the  denominator  of  (67)  is  the  same  as  the  one 

previously  encountered,  and  is  different  from  zero.  Since  the  numerator 

is  of  the  same  form  as  the  denominator;  the  quotient  is  positive.  This 

shows  that  a > 0 when  T > 0 , and  therefore  the  minimum  occurs  when 
o n 

T = 0 s as  was  to  be  proved. 

If  one  of  the  eigenvalues  is  zero,  say  the  last  one,  nothing  is 

changed  in  the  above  work  except  the  equation  for  zn  . Since  the  DE 

for  z is 
n 

z = -w  } (66) 

n n 


the  solution  is 


+ . . o 

z = -w  x + z 
n n n 


It  is  noted  that  if  the  equation  in  the  case  of  distinct  eigenvalues 
is  taken  for  zn  , namely, 


, w X t 

z 7-6  n 

n k 

n 


o + n 
n X_  , 


and  the  limit  is  taken  as  Xn“*-  0 , equation  (89)  is  obtained.  Thus,  it 
is  not  necessary  to  go  through  all  the  algebra  of  successive  substitution. 
All  that  is  necessary  is  to  let  X -^0  everywhere.  This-,  however,  dce3 
not  cause  the  determinants  used  in  the  proofs  of  the  theorems  to  vanish. 
Thus,  the  theorems  are  true  if  one  eigenvalue  is  zero,  and  are  therefore 
true  for  all  cases  where  the  eigenvalues  are  distinct. 

Repeat edEig  envalues 


Consider  the  case  vh er e D , is  of  the  form 


X 1 
X 1 


D 


(90) 
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for  the  present,  it  is  assumed  that  X.  <f  0 . The  P-  and  N-curves  are  given 
by 


+ Dt,  o + s 
z - m ^6  (z  - u)  , 


(91) 


where 


r 
1 1 


^Dt  _ _Xt 


1 t 


o 


t"'1 

• • ; 
*n-2 

U 

(n-2)J 


1 t 


(92) 


The  solution  of  the  minimal  problem  will  follow  the  same  lines  as 
in  the  last  section  and  most  of  the  proof  does  not  need  to  be  repeated. 
In  fact,  all  that  is  necessary  is  to  verify  that  the  two  determinants 
which  correspond  to  (80)  and  (67)  are  different  from  zero. 

Consider  the  determinant  corresponding  to  (87).  Equations  (8ii) 
for  the  path  An  still  hold,  namely, 

, , f -Da  -Da 

zn  1 + (-l)n  \ | 2€  n +...+  (-l)n_x2€-  1 

_T}o  1 

+ (-l)n?“°|u  . (93) 

Differentiating  both  sides  of  this  vector  equation  with  respect  to  a 


gxves 


-Da  -Da  . . -Da. 

0 = -D(  26  nu  -26  n ua  , + o.o  + (-l)n  *€  ~ua, 

n-±  x 

-Da  A 

+ ( -l)n6  0 ua 


(9b) 


T * 

t . 

i. 


0 9 0 


jwy*  ■ < ■ ,r,!fPQfc  W **'  "m  wy  -v* 
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Sir.ca  D is  nonsingular,  the  vector  in  braces  must  be  zero.  Taking  the 
n-compcnent  equations  of  ( 9U)  and  solving  for  , gives,  after  some 
algebra, 


n-l 

r(-*n>k 

Z_y  k:  ak+l 

o 


k l Uk+1 


n-2 

\ '<vk 

kJ  \+2 


(-a1)K 

FI  V-2 


u , - a u 
n-l  n n 


k l k+1  ’ 


u , - a,u 
n-l  1 n 


n-l 

Z_^  ki  \+l 


-Xa 
€ ° 


Fortunately,  these  determinants  can  be  vastly  simplified.  Per- 
forming obvious  elementary  operations  on  rows,  starting  from  the  n 
row  and  working  upwards,  results  in 
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-Xa 


n 


a 

o 

T 


-Xa 


n-1 

n-1 

A A 

a 

a -1 

• 

• 

• 

n 

n-1 

n-2 

n-2 

n-2 

a 

a , 

a-. 

n 

n-1 

1 

• 

m 

• 

• 

• 

• 

• 

• 

• 

a 

n 

a 1 
n-1 

• 

« 

• 

•y 

"1 

1 

2. 

• 

• 

• 

1 

n-1 

n-1 

a 1 

• 

• 

- 

a 

n-1 

0 

• 

• 

• 

• 

• 

• 

• 

• 

a . 

a 

n-1 

0 

1 

• 

• 

« 

• 

1 

(96) 


The  determinants  appearing  are  the  well-known  Vandermonde  deter- 
minants and  are  different  from  zero  if  the  elements  a^  are  distinct, 
which  is  the  caije  here.  Since  the  numerator  and  the  denominator  have 
the  same  form,  the  quotient  is  positive,  as  was  to  be  proved. 

In  the  same  way,  the  determinant  that  corresponds  to  (80)  can  be 
shown  to  be  different  from  zero.  In  fact,  it  turns  out  to  be  the 
Vandermonde  determinant  again. 

Therefore,  the  main  theorem  holds  when  the  matrix  D is  the  ele- 
mentary Jordan  matrix  (90)  and  the  eigenvalue  is  different  from  zero. 

If  the  eigenvalue  in  (90)  is  zero,  then 


(97) 
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The  equations  of  the  P~  and  ;M-curves  can  be  obtained  easily  by 
solving  z - Dz  - w in  the  ordinary  manner.  They  are 

z = -Iw-t  + v,.  —*  ...  + w j+  [z°  + z°  t + ...  + z° 

1 \ 1 c Ji  n ni  / i 1 k n (n-l)il 


, = - ■ w t + w |t  ) ♦ (z°  + z°  t) 

n-1  \ n-1  n 21  j n-1  n 


+ . . o 

z = - w t * z 
n n n 


(98) 


J 


Again,  it  turns  out  that  (98)  can  be  obtained  from  the  solutions 
(91) , where  X / 0,  simply  by  taking  the  limit  as  X 0 . Since  -he 
determinants  in  (96)  do  not  involve  X , the  determinants  obtained  for 
X = 0 will  be  exactly  the  same.  Thus,  the  main  th eorem  also  holds  for 
an  elementary  Jordan  matrix  when  the  eigenvalue  is  zero. 

The  General  Case 

Suppose  the  matrix  D has  the  form 


A 


U = 


lA* . O 


(99) 


where  each  of  the  A ^ is  an  elementary  Jordan  matrix  of  the  form 


j xi  ] 


Ai 


xi 1 o 

9 « 

0 9 

O f 


(ICO) 


The  matrices  can  be  one  by  one  matrices  so  that  purely  diagonal 

boxes  can  be  included  in  D . 
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The  system 

z = Dz  - w (101) 


will  be  nondegenerate  if  w^  / 0 , i = 1,2,..., n and  if  an  eigenvalue 
appears  in  only  one  matrix  of  the  form  (100). 

To  prove  that  the  main  theorem  holds  in  this  general  case,  it  is 
only  necessary,  again,  to  consider  the  determinants  corresponding  to  the 
two  determinants  used  before.  In  order  to  avoid  too  many  rotational  dif- 
ficulties, the  following  special  case  is  considered; 


(102) 


■b 


y (103) 


u II H I rniinn mun-m nr  t r wTrmrwrm “irr H 111'  ' ffnt1T,lWf‘*li“fftT"*l1  Tff1l"lT'1* ITT W fff'trr fl||ggrj*f^l"'TTM*‘*^T* 
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It  is  assuir.ed  that  A ^ is  of  order  k^  , ^ jx  orc*er  > and  A 

of  order  k^  • It  is  also  assumed  that  are  distinct  and  not 

zero  ( this  is  necessary  lor  the  system  to  be  nondegenerate) . 

Let  z^  and  v be  partitioned  as  follows: 


z = i z 


V = I wIj; 


(lOil) 


where  z ^ is  a k^-dimensional,  z^  a k^ -dimensional,  and  zjjj 
k^-dimensional  vector,  and  the  same  for  w^  , w^  , wtxi  * 

The  DE‘s  under  consideration  can  be  written  as 


z.  = A . 2.  - w. 
i i r l 


(i  = I,  II,  III) 


(105) 


These  are  the  same  type  of  matrix  equations  that  have  been  solved  pre- 
viously. After  a little  reflection,  the  relevant  determinant  can  be 
written  without  further  calculation.  The  denominator  of  the  determinant 
used  in  the  proof  of  the  main  theorem  is  (106)  below: 


-X-  a , 
, 1 n-1 


-X  a 

, 1 o 


c'Vn-l 


>„  i Vi 

n-J  • \ k 

v a ,) 
n-x 


-Xv  a 
e kl  o 

-Xa  k_-l 

€ 0 a 2 

o 
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The  only  comment  that  is’  necessary  concerns  the  appearance  of  such 
factors  as  n“^  vhich  did  not  appear  when  the  special  case  of  the 

matrix  Ajj  was  considered  previously.  The  reason  is  that,  previously, 
such  common  factors  of  each  column  cf  the  determinant  could  be  taken  out 3 
this  cannot  be  done  in  the  general  case. 

The  fact  that  this  determinant  is  different  from  zero  follows  from 
problem  75  ana  the  solution  to  problem  76  in  part  V cf  Reference  3.  The 
rest  of  the  arguments  given  previously  apply  without  change,  so  that  the 
main  theorem  remains  true  in  this  case. 

Clearly,  any  number  of  elementary  Jordan  matrices  appearing  in  D 

■JJ- 

will  yield  to  the  same  treatment.  Therefore,  the  main  theorem  holds  true 
for  any  nondegenerate  system! 

Statement  of  Results 

The  original  problem  posed  ha3  now  been  solved  when  the  matrix  of 
the  system  has  real  eigenvalues.  The  result  can  be  stated  quite  simply 
as  follows: 

Theorem  15*  Consider  the  n-dimensional  system 

x = Ax  - e (e  / 0)  , 

where  A has  real  eigenvalues.  If  a_  path  exists  from  any  point  in  the 
phase  space,  then  a unique  minimal  path  exists  and  is  the  path  that  con- 
tains a minimum  number  of  corners. 

The  above  theorem  holds  equally  well  for  degenerate  or  nondegenerate 
systems.  The  only  difference  is  in  the  region  of  existence  of  minimal 
paths. 

For  a nondegenerate  system,  minimal  paths  exist  for  all  points  in 
some  n-dimensicnal  neighborhood  of  the  origin.  If  the  eigenvalues  are  all 
nonpositive,  minimal  paths  exist  from  all  points  in  the  space.  If  the 
eigenvalues  are  all  positive,  minimal  paths  exist  in  a bounded  neighborhood 
of  the  origin.  If  some  of  the  eigenvalues  are  positive  and  some  nonposi- 
tivs,  minimal,  paths  exist  in  an  unbounded  region  lying  between  pairs  of 
parallel  hyperplanes. 


* Theorem  12. 


For  a degenerate  system,  ail  the  minimal  paths  lie  on  some  hyperplane 
through  the  origin.  The  remarks  in  the  preceding  paragraph  hold  for  this 
hyperplane. 
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CON CLUDING  REMARKS 


Some  matters  which  are  pertinent  to  practical  applications  of  the 
results  are  briefly  discussed  belcw: 

if  a control  system  is  to  be  designed  for  a minimum  time  of  cor- 
rr-.cti'.-'  of  an  error  and  its  derivatives,  the  results  in  this  report  and 
those  of  Reference  1 show  that  the  best  control  is  a properly  designed 
limit  control  or  bang-bang  control.  The  mathematical  features  of  the  de- 
sign of  such  a control  have  been  demonstrated  in  many  cases. 

In  the  practical  design  of  such  a system , although  the  control  me- 
chanism itself  is  as  simple  as  possible  (a  switch  or  relay),  a computing 
device  must  be  used  along  with  it  to  instruct  the  control  mechanism  when 
to  switch  from  one  value  to  another.  Before  such  computers  are  designed, 
a great  deal  of  numerical  computation  must  be  done;  however,  the  computa- 
tions are  not  too  complex.  They  would  involve  only  the  tabulation  of  the 
solution  of  certain  linear  differential  equations. 

It  should  be  mentioned  that  a servo-motor  control  system  of  the  type 
discussed  has  already  been  built  (see  Reference  U).  The  differential 
equation  of  the  system  was  one  of  the  simpler  ones,  namely,  x = -1  . The 
system  was  found  to  give  excellent  results  and  to  be  much  lighter  and 
occupy  less  volume  than  a standard  type  of  system  under  the  same  power 
conditions. 

It  should  also  be  mentioned  that  the  design  of  a limit  control  system 
depends  on  the  parameters  of  the  system  to  be  controlled.  However,  it 
certainly  cannot  be  expected  that  a control  can  be  designed  that  will  give 
optimum  performance  without  taking  into  account  the  particular  parameters 

of  the  svs+pm  t.r.  ho 
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